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Abstract 

The Gamma-Gamma (GG) distribution has recently attracted the interest within the research commu- 
nity due to its involvement in various communication systems. In the context of RF wireless communi- 
cations, GG distribution accurately models the power statistics in composite shadowing/fading channels 
as well as in cascade multipath fading channels, while in optical wireless (OW) systems, it describes 
the fluctuations of the irradiance of optical signals distorted by atmospheric turbulence. Although GG 
channel model offers analytical tractability in the analysis of single input single output (SISO) wireless 
systems, difficulties arise when studying multiple input multiple output (MIMO) systems, where the 
distribution of the sum of independent GG variates is required. In this paper, we present a novel simple 
closed-form approximation for the distribution of the sum of independent, but not necessarily identically 
distributed GG variates. It is shown that the probability density function (PDF) of the GG sum can be 
efficiently approximated either by the PDF of a single GG distribution, or by a finite weighted sum of 
PDFs of GG distributions. To reveal the importance of the proposed approximation, the performance of 
RF wireless systems in the presence of composite fading, as well as MIMO OW systems impaired by 
atmospheric turbulence, are investigated. Numerical results and simulations illustrate the accuracy of the 
proposed approach. 
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I. Introduction 

In communication theory, channel statistical modehng is very crucial, since it can be applied in the 
design and performance evaluation of various communication systems. A distribution which has recently 
attracted the interest within the research community due to its involvement in various communication 
systems, is the so-called Gamma-Gamma (GG) distribution. This distribution is equivalent to the squared 
Generalized-K (Kg) distribution [I] and can be derived from the product of two independent Gamma ran- 
dom variables (RVs). Moreover, for certain values of its parameters, it coincides with the if-distribution, 
which in the past has been widely used in a variety of applications, including the statistical characterization 
of the intensity of SAR images [2], as well as to model the statistics of the reverbation intensity in 
underwater communications [3]- [4]. Of particular interest is the application of the GG distribution 
in optical wireless (OW) systems, where transmission of optical signals through the atmosphere is 
involved. In these systems, a major performance hmiting factor is the turbulence induced fading, i.e., 
rapid fluctuations of the irradiance of the propagated optical signals caused by atmospheric turbulence, 
which can be accurately modeled using the statistics of the GG distribution [5]. Furthermore, in recent 
years, GG distribution has also been apphed in the field of RF wireless communications; specifically, 
to model the power statistics in composite fading channels [1], [6]. Additionally, since it includes the 
well known Double-Rayleigh model [7], GG distribution can be further employed to model the power 
statistics in cascade multipath fading channels (which occur, e.g., in keyhole or in mobile-to-mobile 
communication scenarios [8]). The reason for employing GG model in RF wireless systems, is that this 
distribution is general enough to accurately capture the effects of the combined shadowing and multipath 
fading or cascade multipath fading which are encountered in mobile communication channels. 

Although GG channel model is analytically tractable in the performance analysis of various single input 
single output (SISO) wireless communication systems [1], [9]- [12], difficulties arise when studying the 
performance of certain diversity schemes in multiple input multiple output (MIMO) wireless systems. 
Specifically, these difficulties appear when the distribution of the sum of independent GG variates is 
required and have their origin in the fact that a straight derivation of this distribution is analytically 
infeasible, due to the involvement of the modified Bessel function of the second kind. 

In the past, a hmited number of pubhshed works deaUng with the distribution of the sum of independent 
GG variates and its apphcation in communications systems, appeared in the technical literature. In the 
context of mobile communications, where GG distribution models the statistics of the signal-to-noise 
ratio (SNR) in the presence of Kg fading channel model, the maximal ratio combining (MRC) receiver 
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has been investigated in [13] and [14]. However, in [13], the sum of independent GG variates was not 
investigated (shadowing term was common in every diversity branch), while in [14] the expressions that 
were derived for the statistics of the SNR at the output of the combiner, do not holcj]. In the context of 
OW communications, the statistics of the sum of GG variates have been used in the performance analysis 
of MIMO systems operating over GG turbulence model and employing equal gain combining (EGC) at 
the receiver. In [15], the performance of such a system was investigated over identical OW links, using an 
infinite power series representation for the probability density function (PDF) of the turbulence-induced 
fading term at the output of the receiver. Although this approach allowed the derivation of simple and 
accurate expressions at the high SNR regime, it was not computationally attractive when the number of 
the transmit/receive apertures increased and/or the underlying OW links were non identically distributed. 

In this paper, we address to this cumbersome statistical problem by applying a novel and simpler ap- 
proach. We present novel closed-form expressions that approximate efficiently the PDF of the distribution 
of the sum of independent GG variates based on the following well known issues: a) GG distribution is 
derived from the product of two independently distributed Gamma RVs, and b) the distribution of the 
sum of Gamma variates is analytically tractable. Our analysis encompasses the case where the variates 
involved in the sum are identically distributed, as well as the case of non identical variates. Furthermore, 
in order to reveal the importance of the proposed statistical formulation, we study the performance of 
two different systems; an RF wireless system operating over the fading model and employing MRC 
diversity scheme at the receiver and a MIMO OW system operating over strong turbulence channels and 
employing EGC at the receiver. For these systems, closed-form expressions that approximate significant 
system performance metrics, such as bit error rate (BER) and outage probability, are obtained. 

The remainder of the paper is organized as follows. After a brief introduction of the GG distribution 
in Section |II1 novel closed-form expressions that approximate the PDF of the sum of GG variates are 
obtained in Section |llll Moreover, in Section |IV] the obtained results are applied to derive closed-form 
expressions for the performance evaluation of an MRC receiver operating over the fading model, 
while in Section |V] the same results are further applied to the performance evaluation of MIMO OW 
systems operating over strong turbulence channels and employing EGC at the receiver. Finally, in Section 
[Vll useful concluding remarks are provided. 

'in [14, Eq. A-2], it is assumed that the parameter a = fc — m is not an integer. However this comes in contradiction with 
the assumption that the parameters k and m in [14, Eq. A-4] are integers, which results to the incorrect expressions of [14, Eq. 
11] and [14, Eq. 12]. 
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II. The GG Distribution 



Let 7^0, then the PDF of a three-parameter GG RV, which is derived from the square of a Kg- 
distributed RV, is given by [1] 



f-y {i;k,m,n) 



2[Km) 2 -y 2 



k—m 



(1) 



T{m)T{k)n— 

where A; ^ and m ^ are the distribution shaping parameters, (•) is the modified Bessel function 
of order u [16, 8.407/1], r(-) is the Gamma function [16, 8.310/1] and Q, is related with the mean as 
E [7] = $7, with E [•] denoting expectation. 

The distribution in ^ is generic, since it describes various models frequently used in communication 
systems, for several combinations of k and m. Hence, as k —>■ 00, it approximates the well-known Gamma 
distribution (or equivalently squared Nakagami-m [17]), while for m = 1, it coincides with the statistics 
of a squared i^-distributec|^ RV with PDF given by 



f-y (7; ^, ^) 



fc+l k-l 
2k 2 J 2 

fc+1 ^k—m 



k 



(2) 



T{k)n- 

Furthermore for the special case of k = 1 and m = 1, it reduces to the power statistics of the double 
Rayleigh model, frequently used in cascade multipath fading channels, with PDF given by [7]- [8] as 



The n-th moment of 7 is given as [1] 



E[7"] = r 



jF (A; + n) r (m + n) 



(3) 



(4) 



r (k) T (m 

where ^ = Moreover its cumulative density function (CDF) has been expressed using [10, Eq. 7] 
and [16, Eq. 9.31/5] as 

1 



(7; k, m, 0) 



I 1,3 





1 








k, m, 



(5) 



r (k) T (m 

where G [•] is Meijer's G function [16, Eq. 9.301]. 

It is important to note that the GG distribution can be derived from the product of two independent 
RVs, X and y as [18] 

7 = xy (6) 



^In [4] and [5], the statistics of tiie square of a if-distributed RV are also referred as A'-distribution. 
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when are both Gamma distributed with PDF given by 



fi {i; mi,r]i) = — - exp \ i = x,y (7) 

and parameters {mx = k,rix = 1/k) and {niy = m,r]y = ^/m) respective!}!]. 

III. Efficient Approximation to the Sum of GG variates 

Let us consider L independent GG variates denoted by {7;}^^, each having shaping parameters ki 
and mi, and mean The sum of L GG variates, S^, is defined as 

L L 

S-y = "^71 = ^xiyi, (8) 
1=1 1=1 

where xi and yi are Gamma RVs with parameters {ki,l/ki) and {mi,fli/mi) respectively. Eq. dD can 
be rewritten as l l 

(j2 (^E yi^ ^ L~i L 
5*7 = — — Z ^ ^ ~ ^^'^ '■^^ ~ ^^'^ ■ 

i=l j=i+l 

A. Identical Variates 

When the variates of the sum in ([8]l are independent and identically distributed (i.i.d.) (i.e. ki = k, 
nil = rn, Qi = il), {xi}fL^ and {yi}fLi are also identically distributed. Hence, according to the 
unknown distribution of can be approximated by the distribution of the RV S^, which is defined as 



» S, = \ ' . (10) 

with the approximation error, e, given by 

^ L-l L 

{xi-Xj){y,-yj). (11) 

Equivalently, (fTOl l can be written as the product of two RVs s\ and S2, i.e. 

5*7 = S1S2, (12) 

where 

1 ^ 

^It follows from symmetry that it is equivalent to consider the set of parameters (rux = k,rix = ^) and (^niy = m,riy = 
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and 

L 

S2 = ^yi- (14) 

1=1 

Since the sum of i.i.d. Gamma variates remains Gamma distributed [19], it can be proved that si and 
S2 are both Gamma distributed with set of parameters [Lk, and (Lm, ^) respectively. Hence, 
according to S-y will be GG distributed with shaping parameters 

k^_^ = Lk, (15) 

= Lm, (16) 

and mean 

ng= in. (17) 

The accuracy of the proposed approximation depends on the approximation error defined in (fTTI ). The 
exact PDF of the error is difficult to be derived due to its complicate definitiorB However, its first 
moments, which are also indicatives of its statistical behaviour, can be calculated using (ITTT) . Specifically, 
according to the Appendix, the mean of e is equal to 0, while its variance depends from k, m, $7, and 
L, according to 

EH=(L-l)lii. (18) 

It is obvious from the above equation that the variance of the approximation error increases for a certain 
combination of k, m and il, as the number of the RVs of the sum in ([8]l increases, having as a result 
the approximating distribution of Sy to loose its accuracy. 

In order to improve the accuracy of the proposed approximation, an adjustment parameter is introduced 
that modifies the shaping parameters of the approximating distribution of Sy. Specifically, we assume 
that the maximum of the shaping parameter^ k^ is modified by an adjustment parameter, e-y, according 
to 

kg =Lk + e-y. (19) 
The adjustment parameter, e^, is evaluated through the following optimization problem 



arg mm 

£-, 



E 



QU 

Jy 



1,...4 (20) 



''it was observed that it converges to the Laplacian distribution as the number of the RVs of the sum in ([8) increases. 
^Since K-i, (x) = {x), we assume for convenience and without loss of generality that ks. ^ rrig . 
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where E 



are the moments of the distribution of S^, which can be derived from Q using the 
parameters {^kg ^, as provided by ( fT9l ). (fT6l ) and (fTTl ) respectively; E [S*^] are the moments 

of the distribution of S^, calculated using the multinomial expansion, according to 

E [^.1 = E E - E (:) (;;;)-(;;r:)^ [^i""] ^ [^^""^^] [^^^-^ ^^^^ 

!^l=0!/2=0 i^I,_l=0 ^ ^ \ ^ \ ^ 

and using Q for the respective parameters {k, 171,0.). 

The optimization problem in (l20b is a nonlinear multiple function problem, which is difficult to solve 
analytically, yet not impossible to derive an approximative solution numerically. After applying non- 
linear regression methods [20], it was found that the adjustment parameter depends on L, k and m with 
a function of the form of 

, , , -0.127 -0.95A;- 0.0058m 

e^(L,k,m) = (L-l) . (22) 

^ ^ ^ 1 + 0.00124A; + 0.98m 

Hence, using (l22l ) in conjunction with ([T9l ). ([T6l ) and ( fTT] ). the parameters of a single GG distribution are 

defined, which accurately approximates the distribution of the sum of L i.i.d. GG variates. 

B. Non-Identical Variates 

When the GG variates of the sum in ^ are independent and not identically distributed (i.n.i.d.), but 
have one shaping parameter in common, as it happens in most practical applications, (i.e. ki = k, but 
mi and Ui are different), the unknown PDF of the sum can still be approximated by the PDF of the RV 
Sy, as defined by ([T2l ). 

As in the i.i.d. case, Sy can be written as the product of two RVs si and S2, defined by ([T3] ) and ([141 ) 
respectively. Since ki = k, {xi}^^^ are i.i.d. and si is Gamma distributed with parameters (LA;, 
However, the derivation of the distribution of S2 is not straightforward, since {yi}^^ are not identically 
distributed. In order to derive the PDF of S2, the exact closed-form expressions for the sum of non- 
identical Gamma variates presented in [21] are used. According to this approach, the PDF of S2 can be 
written as a nested finite weighted sum of Gamma PDFs, 

L rui / O \ 

(^) = E E 0' > mti) fy Uj, — ), (23) 

i=lj=l ^ 

where y is a Gamma distributed RV with PDF defined by (jV]) and the weights can be easily and quickly 
evaluated using the recursive formula of [21, Eq. (8)] for the parameters {mi}jL^ and according 
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to 



q=l j = l 



(24) 



where t = 1, rui — 1 and with 



llh=i i=i ^ 



3+1- 



The PDF of the product of s\ and is evaluated as 



Using (I23]) and [16, Eq. 3.471/9], Eq. (Eg) yields as 

L m. 



(^) = X^X^^^ (^-J^i^T-J^i >{f^/}f=i) A iz]Lk,j, — j , (27) 

where 7 is a GG distributed RV with PDF defined by ([T]). Hence, an efficient approximation to the PDF 
of the sum of L non-identical GG variates, when one of the shaping parameters remains the same for 



all variatefl can be a nested finite weighted sum of GG PDFs. 

IV. Application in RF Wireless Systems 

A. System Model 

Let us consider a diversity receiver with L branches operating over the composite Kg fading channel 
[14]. The equivalent complex baseband received signal at the /th (/ = 1, 2, L) branch is given by 

zi = shi + ni (28) 

where s is the transmitted complex symbol with energy E and |-| denoting absolute value, hi is 

the channel's complex gain in the path between the transmitter and the Ith branch, and n; is the complex 
Additive White Gaussian Noise (AWGN), having single sided power spectral density No and assumed 
to be identical in all branches. 

*Note that due to symmetry, the same approximation also holds when rrii = m and kt are different, by interchanging ki and 
TUi in l |27| l. i.e. rrii = ki and k — m. 
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Since operation over Kg- fading channel model is considered, the square of the fading envelope, 
Rf = \hi\^, is statistically described by the PDF of Eq. ([T]) with shaping parameters ki and mi, and mean 
E [Rf] = ^i- It follows that the instantaneous SNR of the Zth receiving branch, which is defined as 

-» - % 

is also GG distributed with PDF given by ([Til, shaping parameters equal to ki and mi, and mean equal 
to the average input SNR of the branch defined by 

7, = (30) 

Furthermore, maximum ratio combining (MRC) technique is applied at the receiver and hence the total 
SNR per symbol at the output of the receiver is 

7T = E7. = fE^/- (31) 

1=1 " 1=1 

In the analysis that follows, it is assumed that both shadowing and multipath fading effects are 
independent among the diversity branches, i.e. macrodiversity is studied. As a consequence, the diversity 
branches are considered independent. Moreover, since shadowing occurs in large geographical areas, it 
is further assumed that the parameter that statistically describes the channel's shadowing effects, i.e. 
ki, remains constant among the diversity branches, i.e. ki = k [14]. Note that the assumption of the 
uncorrelated diversity branches can be further applied in the scenario of diversity reception in cascade 
multipath fading channels, when a rich scattering radio environment is considered (see [8] for examples). 

B. Error Analysis 

The average BER of the under consideration RF system can be evaluated directly by averaging the 
conditional BER, (7), which depends from the type of modulation, over the PDF of jt, f-fx (7)' i-^- 

/•OD 

Pbe = / Pe (7) U (7) dj. (32) 

Jo 

Without loss of generality, two types of modulation are considered, binary phase shift keying (BPSK) 
and differential binary phase shift keying (DBPSK). For BPSK and for high values of the average input 
SNR, the conditional BER is given by (7) = ^erfc {^/^), where erfc(-) is the complimentary error 
function [16, Eq. (8.250/1)], while for DBPSK is given by Pe (7) = ^ exp (—7) [22]. Hence, the average 
BER of the diversity system can be evaluated by 

— I r"^ 

Pke = -J^ erfc ( V7) U (7) d-f (33) 
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Pbe = l I exp (-7) /^^ (7) dj (34) 



and 

Pk„ = 

2 JO 

when BPSK and DBPSK modulation schemes are employed respectively. 

1) Independent and Identically Distributed Diversity Branches: When the received signals at the 
diversity branches are independent and identically distributed (i.i.d.), i.e. nii = m and 7^ = 7, the 
approximation for the PDF of the sum of i.i.d. GG variates can be used in order to evaluate the integrals 
in (I33] ) and (l34l ). Specifically, the PDF of can be approximated by the PDF of a single GG variate 
with parameters defined by ( fT9l ). (fT6l) and (fTTl) . i.e. 

fiT (7) ~ /t (7; ^T, rriT^^T) (35) 

where kx = Lk + e^^ , uit = Lm, = Lj and e^^ is defined from ((22l) for the set of parameters 
of {L,k,m). By substituting (1351 ) in (l33l ) and (l34b . and using the corresponding BER expressions of 
a Single Input Single Output (SISO) system with parameters (fcy, m^, 7ji), analytical expressions that 
approximate the average BER of the diversity system can be derived. Hence, using [1, Eq. (8)], the 
average BER for BPSK modulation can be approximated by 

-Gi 



2^T [kr) r (mr) 



l-/3r 

2 ' 2 



2 ' 2 ' 2 



(36) 



where = (3t = kx + rriT — 1 and ax = kr — niT- In a similar manner, using [1, Eq. (9)], the 

BER performance of the DBPSK diversity system is evaluated as 

_ I St /£^\ 

Pbe ~ TT^T exp ^ ^ (Ct) , (37) 



2'^ V 2 , 

where Wx,f, (•) is the Whittaker function [16, Eq. 9.220]. 

2) Independent, but not Necessarily Identically Distributed Diversity Branches: When the received 
signals at the diversity branches are independent, but not necessarily identically distributed (i.n.i.d.), i.e. 
mi and 7^ are different among the diversity branches, the approximation for the PDF of the sum of i.n.i.d. 
GG variates can be used in order to approximate the average BER of the under consideration diversity 
system. Hence, according to (|27] ). the PDF of 7t can be approximated by a nested finite weighted sum 
of GG PDFs, i.e. 

L rrii , . \ 

U (7) ^ ^ H (i, j) A ( 7; ^A;, i, ) , (38) 
,-1 ,-1 V 



i=\ j=l 



where 



=WL {i,j,{m}f=i,{li}ti) (39) 
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are evaluated using (I24b and (I25l ). By substituting the above sum in (1331 ) and (I34b . the average BER of 
the system for BPSK modulation is approximated by 



L nii 



1 



i=l j=l 



while for DPSK modulation, by 

L rrii 
i=l j=l 



exp(-7)/^ r,Lk,j 



d'y. 



(40) 



(41) 



The integrals in the above equations can be evaluated using the corresponding BER expressions of a 
SISO system with parameters (^Lk,j,^^. Hence, using [1, Eq. (8)], the analytical expression that 
approximates the average BER of the BPSK diversity system yields as 

i-ft ft 



L nii 



rpt2\/^r(LA:)r(j) 



G 



2,2 
2,3 



2 ' 2 

.ft+1 



(42) 



where 



Lknii 



^\ and /3j = LA; + j — 1. Moreover, using [1, Eq. (9)], the average BER of 



the DBPSK diversity system is derived as 



(43) 

i=l j=l 

C. Outage Probability 

Outage probability is defined as the probability that the output SNR of the under consideration diversity 
system falls below a specified threshold 7^^ , which represents a protection value of the SNR above which 
the quality of the channel is satisfactory. Hence, outage probability can be evaluated by 



7th 



/7t (7) d'y. 



(44) 



Pout = Pr (7T <lth) = 

1) Independent and Identically Distributed Diversity Branches: When the received signals at the 
diversity branches are i.i.d., i.e. mi = m and 7^ = 7, the PDF of 77' can be approximated by a single 
GG distributed variate with parameters (/cy, mj-, 7^.) according to (1351 ). Hence, the outage probability of 
the diversity system can be approximated by 



7tfc 

Pout ~ / f-y (7; kr, mT,^T) dj, 
lo 



(45) 



which is equivalent to the outage probability of a SISO system with parameters (A;^, mr, 7t)- Using 
(l45l) yields 

1 



Pr 



out 



V{kT)T{mT) 



^1,3 





1 


iTith 






kTjinj^, 



(46) 



where = 

7t 
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2) Independent but not Necessarily Identically Distributed Diversity Branches: When the received 
signals at the diversity branches are i.n.i.d., i.e. mi and 7^ are different among the diversity branches, 
the PDF of 7J' is approximated by a nested finite weighted sum of GG PDFs, according to ( [38l ). Hence, 
the outage probabiUty can be approximated by 



which is equivalent to a nested finite weighted sum of outage probabilities of SISO systems, each having 
the parameters (Lk,j, Using ©, m yields 



where ii = idMBu._ 

D. Numerical Results and Discussion 

Figs. [T] and |2] illustrate the average BER and outage probability of MRC receivers operating over the 
Kg fading model, when the diversity branches have the same shaping parameters and average branch 
input SNR, i.e. ki = k, mi = m and 7; = 7. Approximative analytical results, using Eq. (l36l ) and (|37] ) 
for the BER evaluation and ( |46l ) for the outage probability evaluation, are plotted in comparison with 
Monte-Carlo (MC) simulation results, for an arbitrary number of diversity branches and assuming a certain 
combination of shaping parameters (k = 2 and m = 5). It is observed that there is an excellent match 
between simulation and the approximative results for every input SNR and normalized outage, ^thM, in 
both performance metrics. It is also clearly depicted that the approximative analytical expressions remain 
accurate, even when the number of the diversity branches increases. 

In Figs. |3]|6j the BER and outage probability performance metrics of MRC receivers are depicted, when 
the diversity branches are i.n.i.d., i.e. mi and 7^ are different among the diversity branches. Specifically, 
it is assumed that the average input SNR of Ith branch is given by 7^ = 7^^ exp [—5 {I — 1)], where 71 
is the average input SNR of the first branch and 5 is a decaying factor. Using ( |40l ) for BPSK and (|4T]) 
for DBPSK modulation. Figs. [3l|4] present the approximative analytical results for the average BER, as 
a function of the first brach average input SNR, for several combinations of shaping parameters k, mi 
and decaying factors 5. It is obvious from the figures that the approximative analytical expressions for 
the average BER are close to simulation results (their difference is not greater than 3 dB at target BERs 
equal to 10^^), even when 5 increases. Moreover it is evident that the proposed approximation acts as a 




L m.. 



(47) 





(48) 
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lower bound for high values of SNR and the smaller the number of diversity branches and/or 6 are, the 
more accurate is the bound. Similar behaviour has been also observed for the outage probability, which 
is depicted in Figs. |5]and|6]as a function of the first branch normalized outage threshold, 7th/7i> for 
the same combinations of k, mi and 6. The difference between analytical results, derived from (l48l) . and 
simulation results also lies within 3dB in all cases examined, and the proposed approximation acts as a 
tight lower bound. 



A. System Model 

Consider a Multiple Input Multiple Output (MIMO) optical wireless (OW) system where the infor- 
mation signal is transmitted via M apertures and received by N apertures over strong atmospheric 
turbulence conditions. For the OW system under consideration, it is assumed that the information bits are 
modulated using On-Off keying (OOK) and transmitted through the M apertures using repetition coding 
[23]. Moreover, a large field of view is considered for each receiver indicating that multiple transmitters 
are simultaneously observed by each receiver. This actually leads to the collection of larger amount of 
background radiation which justifies the use of the AWGN model as a good approximation of the Poisson 
photon counting detection model [24]. Hence, the received signal at the qth receive aperture is given by 



where x G {0, 1} represents the information bits, t] is the optical-to-electrical conversion coefficient and 
Vq is the AWGN with zero mean and variance = No/2. 

The term Ipg denotes the fading coefficient that models the atmospheric turbulence through the optical 
channel between the pth transmit and the qth receive aperture. Since operation under strong atmospheric 
turbulence conditions is assumed, according to [9], the parameter which represents the effective number 
of small scale scatterers can be considered equal to 1. Hence, the optical channel in the p-qth transmit- 
receive pair can be statistically described by a GO distribution with parameters k = 1, m = Upg and 

= E [Ipg] [9], where Upg is related to the effective number of large scale scatterers. Furthermore, 
it is assumed that the statistics of the fading coefficients of the underlying channels are statistically 
independent; an assumption which is realistic by placing the transmitter and the receiver apertures just a 
few centimeters apart [25]. 

At the receiver side, the received optical signals from the A'^ apertures are combined using equal gain 
combining (EGC), which is an efficient combing scheme in OW systems [24]- [25]. Hence, the output 



V. Application in Optical Wireless Systems 



M 




(49) 



p=i 
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of the receiver is 

N N M 

7=1 q=l p=l 



Note that a scaling factor of MN appears in ( pOl) . The factor M is included in order to ensure that the 
total transmit power is the same with that of a system with no transmit diversity, while the factor 
ensures that the sum of the N receive aperture areas is the same with the aperture area of a system with 
no receive diversity. 

The received electrical SNR of the OW link between the p transmit and q receive aperture, can be 
defined as [26] 

hpn = V^, (51) 



while its average as fipg = ^ -^^p"! . According to the above definitions, the electrical SNR of the 
combined signal at the output of the receiver, becomes 

= apnWo' ^^^^ 

where h = Ef=i Ep=i W 
B. Error Analysis 

The BER probability of the MIMO OW system under consideration, assuming perfect Channel State 
Information (CSI), is given by [24] as 



N M 



q=l p=l 



where fi{I) is the joint PDF of the vector / = {In, I12, ■■■Imn) of length MN. Furthermore, Q {■) 
is the Gaussian-Q function defined as Q (y) = (l/\/27r) JJ^ exp (— dt and related to erfc(-) by 
erfc (x) = 2Q (\/2x) . Equivalently, Eq. ([53] ) can be evaluated as 

Pe = l [ fi^ (!) erfc ( ^ 1 I dl (54) 

2 70 \2V2NMa^ J 

where //^ (/) is the PDF of It- 

1 ) Independent and Identically Distributed OW Links: When the turbulence induced fading coefficients 
of the underlying optical links of the MIMO system are independent and identically distributed, i.e. 
Qpq = a and E [Ipg] = lo, the PDF of It can be approximated by the PDF of a single GG variate, i.e. 

fiAl)- fj{I;kT,mT,nT) (55) 
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where kx = MNa + e^^ , rriT = MN, Q,t = MNIq and e^^ is calculated from (1221 ) for the parameters 
{MN,a, 1). Hence, the BER probability of (l54l ) is approximated by 



I di. 



(56) 



.2y2NMo-„ 

The integral of ( [56l ) can be solved using Meijer's G-functions and their properties. Hence, by substi- 
tuting the PDF of the GG distribution according to expressing the (•) and the erfc (•) integrands 
in terms of Meijer's G-function according to [27, Eq. (8.4.23.1)] and [27, Eq. (8.4.14.2)] respectively, 
and using [27, Eq. (2.24.1.1)], the BER is expressed as 



V^r{kT)r{mT) ^'^ 



1— fcr 2— fcj l—m-r 2— ■ 
2 ' 2 ' 2 ' 2 ' 







' 2 



(57) 



where jj, denotes the average electrical SNR of each OW link. 

2 ) Independent and Not Identically Distributed OW Links: When the turbulence induced fading co- 
efficients of the underlying optical links of the MIMO OW system are independent, but not identically 
distributed, the PDF of It can be approximated by a nested finite weighted sum of GG PDFs, according 
to ([27]), i.e. 

L rrii / .„ 



1 (58) 

i=i j=i ^ 

where L = MN is the number of the underlying OW links, mi = Upq, ki = 1 and = E [Ipq], when 
p = l, ...M, q = l, ...N and / = 1, ...MN. Furthermore, 



(59) 



are evaluated using (|24|) and (1251) . 

By substituting (1581) to (l54l ). the BER probability of the MIMO OW system is approximated by 

L rrii 



i=i j=i 



rrii 



erfc 



r] 



I dl 



(60) 



The integral in the above equation can be evaluated by expressing its integrands in terms of Meijer's 
G-function, as in the i.i.d. case. Hence, the probability of error is given by 



L rrii 



P. 



where Hi is the average electrical SNR of the ith OW Unk. 



Mi 



1-L 2-L 1-j 2-j 
2 ' 2 ' 2 ' 2 ' 



2 



(61) 
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C. Outage Probability 

As in wireless RF systems, the outage probability is defined as the probability that the SNR of the 
combined signal at the output of the receiver, falls below a specified threshold hth- It is considered as an 
important parameter for OW links to be operated as a part of a data network and is critical in the design 
of both transport and network layer. 

Hence, according to (l52l) . the resulting outage probability of the system is given by 



Pout = Pr (hr < hth) = Pr (It < Ith) 



(62) 



^^Vhth^o is the normalized threshold. 



where Ith 

1 ) Independent and Identically Distributed OW Links: When the underlying channels of the MIMO 
OW system are independent and identically distributed, the PDF of It can be approximated by Eq. (1551) . 
The outage probability of the under consideration system can be approximated by 



Pout^ / fj{I;kT,mT,^T)<iI 
Jo 



(63) 



which is equivalent to the outage probability of a SISO system operating over the GG turbulence model 
with parameters (/ct, "ly, Hence, using Q, a closed-form expression for the outage probability 
yields as 



Pr 



1 



out 



r (kx) r {rriT 



r^2,i 

"^1,3 



kTrriT 



th 



1 



kx, nij', 



(64) 



2 ) Independent and Not Identically Distributed OW Links: When the underlying channels of the MIMO 
OW system are independent, but not identically distributed, the PDF of It can be approximated by Eq. 
(l58l) . Hence, the outage probability of the under consideration system can be approximated by 



A P,L,j 



mi 



dl. 



(65) 



L rrii 

Pout 

i=i j=i 

From the above equation, it is evident that the outage probability of the MIMO OW system can be 
approximated by a finite nested weighted sum of outage probabilities of SISO OW links operating over 



the GG turbulence model with parameters I L, j, 



rrii 



P 



""*^£rfer(L)r(,)'^i'=' 



Using ([S]), an analytical expression is derived as 
Lm,- 



-I 



th 



1 



(66) 
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D. Numerical Results 

In Figs. |7] and [8j the BER and outage probability of MIMO OW systems operating over identically 
distributed strong turbulence channels with parameters a = 4 or a = 10 and Iq = 1, are depicted. 
Analytical results, using (l57l) and (l64l) . are illustrated in comparison with MC simulations. It is observed 
that there is an excellent match between the approximation and the simulations in every SNR regime 
for both performance metrics. It is also clearly depicted that the derived approximative expressions are 
accurate for every MIMO deployment investigated, irrespective the number of transmit and/or receive 
apertures. 

Figs. |9] and [10] depict the BER and outage probability of various MIMO deployments of transmit 
and receive apertures over i.n.i.d. strong turbulence channels, i.e. the underlying OW links have different 
turbulence parameters and different average electrical SNRs. As it is clearly illustrated, the approximative 
analytical results for both performance metrics, derived from (|6TI ) and (l66l ) respectively, are very close 
to the MC simulation results. Specifically, for the MIMO deployments investigated and for practical 
values of average BER and outage probability, the difference between analytical and simulation results is 
not greater than 2 dB. Moreover, it is observed that the proposed approximation acts as a lower bound, 
which becomes less accurate as the number of the underlying i.n.i.d. OW links increases. However, taking 
into consideration that the BER performance metric is difficult or even impossible to be evaluated with 
numerical techniques as the number of the OW links increases [26], the derived closed-form expressions 
can be considered as reliable alternatives to time consuming MC simulations. 

VI. Conclusions 

We examined the statistics of the sum of independent and not necessarily identical GG RVs. Novel 
closed-form expressions were derived that approximated its PDF either with the PDF of a single GG 
distribution, when all the variates of the sum were identically distributed, or with a finite weighted sum 
of PDFs of GG distributions, when the variates of the sum were non identically distributed. Based on the 
obtained statistical formulas, the performance of MRC diversity receivers operating over the Kg fading 
channel, as well as MIMO OW systems operating over strong turbulence channels and employing EGC at 
the receiver, was investigated and major performance metrics were analytically evaluated. The comparison 
between approximative analytical results and simulations demonstrated that the proposed approximation 
is accurate when the diversity branches or underlying OW links are identical, while it serves as a tight 
lower bound for non identical diversity branches or OW links. 
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Appendix 

Since both {a^ijjz^' and {yiljzf are identically distributed squared-Nakagami variates with parameters 
{k, l/k) and (m, J7/m), respectively, their first moments and variances are given by 



E [xi] = 1 
E [yi] = n 

E[yf]-E[mf 



2 1 



m 



(67) 
(68) 
(69) 

(70) 



where I = 1, ...L. 



Due to their independency, the first moment of the the error e can be easily calculated, according to 

L-l L 

E = I E E [x.] - E [x,]) (E [y,] - E [y,]) = 0. (71) 

i=l j=i+l 

Furthermore, the variance of e is equal to its second moment and it is derived from 

2n 



E [a 



1 

l2 



E 



' L-l L 

E E i^i - ^j) iVi - Vj) 

1=1 j=i+l 



(72) 



After expanding the sums in (1721 ) and taking their square, the terms that appear are calculated using (1671 ). 
K (l69l ) and (iTOl ). Specifically, the following terms appear. 



E [{x^ - x]) {xl - xl) {yf - y]) {yl - y^)] 



4£ if i = /land J =5 
if i 7^ /i and j ^ g 
m Hi = hand J ^g 



(73) 



where i,j,h,g = I...L. After some algebra and using (1731 ). (1721 ) is simplified to ([TSl) . 
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Branch Average Input SNR y, dB 



Fig. 1. Comparison of approximate average BER and MC simulation results of MRC receivers with i.i.ddiversity branches 
(k = 2 and m = 5). 
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Normalized Outage Threshold, dB 

Fig. 2. Comparison of approximate Outage probability and MC simulation results of MRC receivers with i.i.d. diversity 
branches (k = 2 and m = 5). 
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Fig. 3. Comparison of approximate average BER and MC simulation results of MRC receivers with i.n.i.d.diversity branches 
(S = 0.5). 
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Fig. 4. Comparison of approximate average BER and MC simulation results of MRC receivers with i.n.i.d.diversity branches 
(5 = 1). 
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Fig. 5. Comparison of approximate Outage probability and MC simulation results of MRC receivers with i.n.i.d. diversity 
branches (5 = 0.5). 
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Fig. 6. Comparison of approximate Outage probability and MC simulation results of MRC receivers with i.n.i.d. diversity 
branches (5 = 1). 
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Fig. 7. Comparison of approximate average BER and MC simulation results for various MIMO OW systems over i.i.d. Unks. 
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Normalized Threshold l_, dB 



Fig. 8. Comparison of approximate outage probability and MC simulation results for various MIMO OW systems over i.i.d. 
links. 
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Fig. 9. Comparison of approximate BER performance and MC simulation results of MIMO systems with M = 2 and N = 1, 
M = 3 and N = 1, and M = 2 and N = 2 transmit and receive apertures. 
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Fig. 10. Comparison of approximate outage probability and MC simulation results of MIMO systems with M = 2 and N = 1, 
M = 3 and N = 1, and M = 2 and N = 2 transmit and receive apertures. 
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